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ABSTRACT 

We have compared the radio emission from a sample of parsec-scale AGN jets as measured by the VLBA at 
15 GHz, with their associated 7-ray properties that are reported in the Fermi LAT 3-month bright source list. 
We find in our radio-selected sample that the 7-ray photon flux correlates well with the quasi-simultaneously 
measured compact radio flux density. The LAT-detected jets in our radio-selected complete sample generally 
have higher compact radio flux densities, and their parsec-scale cores are brighter (i.e., have higher brightness 
temperature) than the jets in the LAT non-detected objects. This suggests that the jets of bright 7-ray AGN 
have preferentially higher Doppler-boosting factors. In addition, AGN jets tend to be found in a more active 
radio state within several months from LAT-detection of their strong 7-ray emission. This result becomes more 
pronounced for confirmed 7-ray flaring sources. We identify the parsec-scale radio core as a likely location for 
both the 7-ray and radio flares, which appear within typical timescales of up to a few months of each other. 
Subject headings: galaxies: active — galaxies: jets — radio continuum: galaxies 



1. INTRODUCTION 

A number of authors have suggested a close connection 
between AGN 7-ray and radio emission, with the strong 
energetics of the jetted relativistic outflow being responsi- 
ble for accelerating particles to the high ener gies needed 
for 7-ray p roduct ion (e.g., Dermer & Schlick eised 119931 : 
ISikoraet alJ 119941: iBeall & Bednareb il999.) . Most of the 
high energy 7-ray sources detected by the EGRET telescope 
of th e Compton Gamma Ray Observatory ( Hartman et al] 
T999h were identified w ith blazars (e.g.. iMattox et al.. ,20011 : 



Sowards-Emmerd et al.l l2003, 2004) which indicates that rel- 



ativistic jets with strong Doppler boosting are dominant sites 
of extragalactic 7-ray production. The EGRET results, how- 
ever, were limited by the sensitivity of the telescope and the 
eiTatic sampling of observations. This did not allow a full 
systematic comparison of the 7-ray and radio properties for a 
complete sample of blazars. 

' Max-Planck-Institut filr Radioastronomie, Auf dem Hiigel 69, 53121 
Bonn, Germany; 

- Astro Space Center of Lebedev Physical Institute, Profsoyuznaya 84/32, 
1 17997 Moscow, Russia 

^ Department of Astronomy, University of Michigan, 817 Denison Build- 
ing, Ann Arbor, Ml 48109-1042, USA; 

* Department of Physics and Astronomy, Denison University, Granville, 
OH 43023, USA; 

^ Dr. Remeis-Stemwarte Bamberg, Universitat Erlangen-Niirnberg, Stem- 
wartstrasse 7, 96049 Bamberg, Germany 

^ Erlangen Centre for Astroparticle Physics, Erwin-Rommel Str. 1, 91058 
Erlangen, Germany 

' CRESST/NASA Goddard Space Flight Center, Greenbelt, MD 20771, 
USA 

^ Universities Space Research Association, 10211 Wincopin Circle, Suite 
500 Columbia, MD 21044, USA 

' National Radio Astronomy Observatory, 520 Edgemont Road, Char- 
lottesville, VA 22903-2475, USA; 

Department of Physics, Purdue University, 525 Northwestern Avenue, 
West Lafayette, IN 47907, USA; 

" Pulkovo Observatory, Pulkovskoe Chaussee 65/1, 196140 St. Peters- 
burg, Russia 

Crimean Astrophysical Observatory, 98409 Nauchny, Crimea, Uki'aine 
" Departament d'Astronomia i Astrofi'sica, Universitat de Valencia, E- 
46100 Burjassot, Valencia, Spain 



iJorstad et al] (12001b carried out 43 GHz Very Long Base- 
line Array (VLBA) monitoring of a sample of EGRET blazars 
and concluded that 7-ray events originate in parsec-scale ra- 
dio knots that move down the jet at apparent superlumi- 
nal speeds. Lahteenmaki & Valtaoia (2003) supported this 
finding on the basis of 37 GHz total flux density single- 
dish monitoring observations of EGRET blazars. They con- 
cluded that the highest levels of 7-ray emission are ob- 
served during the initial or peak stages of radio/mm-wave 
flares. VLBA observations of a radio-sel ected sample of 
bright extragalactic sour ces at 15 GH z (Keller mann et alj 
12004: Kovalev et al. 2005,:iLister & Homan 2 005) have shown 
that jets of EGRET-detected blazars are more compact, con- 
tain faster-moving knots, and are more luminous and more 
linearly polarized than non-EGRET blaz ars. Recent analy- 
sis of a large 6 cm VLBA survey (VIPS, iTavloret al.l 120071) 
has shown that EGRET-detected blazars tend to have higher 
brightness temperatures and greater core dominance. How- 
ever, a direct correlation between radio and 7-ray flu x den- 
sity and luminosity (e.g., |Aner et al. 1996; Miicke et al .ll 19971: 
iLahteenmaki & Valtaoian2003i: iTavlor et all 120071) was not 
found. 

The Fermi Gamma-Ray Space Telescope (previously 
known as GLAST ) was successfully launched in June 2008. 
One of the two instrum ents on board the spac e craft, the Large 
Area Telescope (LAT, lAtwood et al.l 120091) . is an imaging, 
wide-field telescope, covering the energy range from about 
20 MeV to more than 300 GeV Its characteristics allow for 
a quasi-continuous systematic study of the whole 7-ray sky 
with unprecedented sensitivity. The Fermi LAT team has re- 
cently released results of the first three months of observa- 
tions (August - October 2008) in the form of a bright source 
list'"* ( Abd o et alJ2009 b). In this Letter we compare the 7-ray 
and parsec-scale radio properties of the initial (> lOti) LAT 
AGN detections that are positi onally associated with bright 
radio-loud blazars (Abdo et alj,2009b,a ). We call this sam- 
ple throughout the paper 'ELBS' — Extragalactic LAT Bright 
Source sample. 

Throughout this letter, we use the term "core" as the appar- 

http://fermi.gsfc.nasa.gov/ssc/data/access/lat/bright_src_list/ 
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ent origin of AGN jets that commonly appea rs as the bright- 
est feature in V LBI images of blazars (e.g.. lLobanovl[T998t 
lMarscheij|2008l) . We use the ACDM cosmol ogy with Hp = 
71 km s"^ Mpc~\ fim = 0.27, and flA = 0.73 (iKomatsu et al.1 
l2008h . 

2. THE RADIO DATA AND SOURCE SAMPLES 

The MOJAVE sample of 135 extragalactic radio sources 
was rigorously selected on the basis of the parsec-scale flux 
density at 15 GHz (5vlba > 1-5 Jy), and is part of an on- 
going monitoringjgrogram with the VLBA at this frequency 
(for details, see iLister & Hom an 2005; ILister et al.ll2009ah . 
The sample consists of the brightest extragalactic jets north of 
declination -20°, and is dominated by blazars (101 quasars 
and 22 BL Lac objects). Thirty one AGN from this sample 
are reported by the Fermi LA T team in the 3-m onth Fermi 
LAT Bright 7-ray Source List (lAbdo et al.l l2009b). The MO- 
JAVE program also began, in late 2008, to monitor weaker 
LAT-detected AGN (flux density S' > 0.2 Jy with S > -30°), 
some of these were used for the analysis in § 13.11 where 
we supplement the MOJAVE complete sample with 46 other 
LAT-detected AGN with 5' > 0.2 Jy. We emphasize that the 
sources plotted in Figure [T] do not represent a complete sam- 
ple selected on either parsec-scale radio or gamma-ray flux, 
but this Figure does show the best available comparison at 
this time. We a lso note that many LAT-detected BL Lacs 
dAbdo et al.l2009 a) are radio weak and do not enter our radio- 
selected sample and analysis. 

We have attempted to use radio data that are contempora- 
neous with the LAT measurements, to avoid, as much as pos- 
sible, effects due to variability. We use VLBA observations 
made during May - December 2008 for most of the sources. 
For the remaining objects not observed then at the VLBA, 
we use recent data from RATAN-600 observations in Septem- 
ber/October 2008 (e.g., Kovalev et al. 1999) or the U MRAO 
obser vations in August 2008 - January 2009 (e.g., All er et al.l 
l2003h . We infer the 15 GHz parsec-scale flux density with 
the assumption that the total single-dish flux well represents 
the parsec-scale emission. We used previ ous RATAN results 
together with the VLBA calibrator survey jPetrov et al.ll2008l 
and references therein) to confirm this assumption, follow- 
ing the compactness analysis bv Kovalev et al. (200 5). Sev- 
eral sources with no VLBA flux density measurements but 
with strong kiloparsec-scale emission were not included in the 
analysis. If a source was observed more then o nce by RATAN 
or UMRAO, an average value is used. Note that 'K ovalev et al.l 
(|2005) have compared flux density scales of the three pro- 
grams run at the VLBA, the UMRAO, and RATAN, and found 
that they agree within the errors, that typically are 5 per cent. 

In addition to these quasi-simultaneous measurements of 
the radio flux densities, we have used in the analysis in § |3] 
(i) the results from our earlier 15 GHz flux density measure- 
ments within the 2cm VLBA/MOJAVE program presented by 
lOster et al.. (2009a,), and (ii) brightness temperature me asure- 
ments of the jet co res made since 1999 (Kovalev et al .l 120051; 
iHoman et al.0OO6l) and recently updated by McCormick et al.l 
(120091) . We did not restrict the analysis in galactic latitude, b; 
but galactic plane emission does decrease the LAT detection 
sensitivity for \b\ < 10° (Abdo et al. 2009b), and this could 
bias our conclusions. To test this we re-did our analysis for 
1^1 > 10°, and achieved the same results at approximately the 
same confidence level in every single test. We also note that 
the two brightest radio lenses on the sky which were detected 
by Fermi LAT, B 0218+357 and PKS 1830-211, are not in- 
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Fig 1 . — Average Fermi LAT 100 MeV - 1 GeV photon flux lAbdo et all 
l2009b*) versus quasi-simultaneous 15 GHz flux density. Filled circles repre- 
sent total VLBI flux density while open ones — single-dish flux density. The 
single dish flux densities are representative of the parsec-scale emission in 
these objects as described in §|2] 



eluded in the MOJAVE sample (ILister et al.1 l2009al) and are 
not part of the analysis in this paper. 

3. RESULTS 
3.1. Quasi-simultaneous "f-ray vs. radio emission 

Figu re[T]plots the average photon flux measured by the LAT 
dAbdo et al. 2009b) in the energy range 100 MeV - 1 GeV 
during three months in August-October, 2008, against the 
parsec-scale radio flux density, Sq^, measured by the VLBA 
at 15 GHz within several months of the LAT observations, 
namely between May 2008 and December 2008. This is cal- 
culated as the total flux density over the entire VLBA image. 
VLBA data points for 38 sources are supplemented by single- 
dish measurements for 39 objects as discussed in § |2l The 
points in Figure [T]represent the AGN in the ELBS list (> lOcr 
LAT detections) with radio flux density greater than 0.2 Jy and 
declination 6 > -30°. This sample is about two times larger 
than the ELBS-MOJAVE sample. We note that the redshift 
distribution of the latter is not statistically di fferent from that 
of the full ELBS sample ( Lister et al.ll2009bl) . 

The non-parametric Kendall's tau test confirms a positive 
correlation in Figure [T] at a confidence level > 99.9 per cent 
for the 100 MeV - 1 GeV energy band. The test provides 
the same result when restricted to only those sources mea- 
sured by the VLBA (confidence 99.5 per cent) or when re- 
stricted to the radio brightest MOJAVE-ELBS sources (confi- 
dence 99. 1 per cent) or when restricted to ELBS quasars only 
(confidence 9 8.9 per cent). See dis cussion on the possible bias 
by BL Lacs in lAbdo et al] (l2009al) . The same analysis for the 
second LAT energy band, 1-100 GeV, also shows a positive 
correlation, but at a lower confidence level, 86 per cent. 

These results indicate a direct relation between the 7- 
ray and pa r sec-sca le synchrotron radiation. We note that 
lAbdo et aP (l2009al) found a much less significant correla- 
tion using non-simultaneous total radio (VLA) flux density 
in their 7-ray selected sample. In terms of previous EGRET 
results, it is possible that a positive correlation was not found 
by, e.g., .Aller et al.. (1996) ; .Liihteenmaki & Valtaoia (2003.) ; 
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Fig. 2. — Radio 15 GHz flux density distributions for the complete MO- 
JAVE sample. Left panel: mean flux density, Sgg, measured between May and 
December 2008. Right panel: mean flux density, S!,^ Qy, measured between 
January 1999 and September 2007. The shaded areas show LAT-detected 
(ELBS) sources. 



iTaylor et"an (120071) in the EGRET vs. radio analyses because 
of the lack of homogeneous statistics, limited EGRET sen- 
sitivity and/or use of non-simultaneous measurements in the 
two bands. A more thorough analysis with direct implications 
for the emission models will be possible when the LAT 7-ray 
energy flux data become available. 

3.2. 75 GHz parsec-scale flux density and its variability 



The correlation reported in § 13.11 is further supported by 
the comparison of the average radio flux density of the ob- 
jects detected by Fermi to that of the non-detected sources 
for the complete MOJAVE sample. We assume that LAT non- 
detected objects have statistically lower photon flux and there- 
fore should have lower radio flux. We analyzed this with 1 14 
VLBA data points and 21 single-dish measurements made 
in 2008 for the complete MOJAVE sample. The 5J)g radio 
flux densities of parsec-scale jets detected by Fermi (ELBS 
sources) are higher (mean value 3.0 Jy vs. 1.6 Jy) and have a 
different distribution from the non-detected ones as indicated 
by the Kolmogorov-Smirnov (K-S) test at a confidence level 
99.9 per cent (Figure |2]i. 

The MOJAVE database'^ contains a complete set of all cal- 
ibrated VLBA observations of the MOJAVE sources since 
1999 and includes data which we reduced from the pub- 
lic NRAO archive in addition to our own observations 
jLister et al.ll200 9a^. This provides us with an estimate of the 
mean flux density state of each source, S\c)_q-i, over the previ- 
ous eight years, between January 1999 and September 2007. 
We have made the same statistical test as above, flux of LAT- 
detected vs. LAT non-detected, using 599.07. The general con- 
clusion remains the same as indicated by the K-S confidence 
99.5 per cent, although the difference between ELBS and non- 
ELBS sources becomes less pronounced. The difference is 
also found at the 99.6 per cent K-S confidence for distribu- 
tions of the 2008 radio luminosity, L^g, if compared for MO- 
JAVE quasars only (Figure|3]l. BL La cs in the MOJAVE sam- 
ple have statistically lower redshifts jLister et al.ll2009ab and 
were excluded from the luminosity analysis. 

Finally, an intriguing result emerges if we examine the ra- 
dio activity level of the ELBS sources by comparing the flux 
density measurements in 2008, Sq,^, with the respective aver- 
ages over previous eight years, Sgi^-Q-i- We define the 2008 
parsec-scale activity index as Vq^ = (Sq^ - S\q_Q-i) / Sqq_Q-,. The 
ELBS list includes a flag that indicates whether or not the 
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Fig. 3. — Distribution of 15 GHz luminosity calculated from SJjg (Figure|4) 
for the quasars in the MOJAVE sample. The shaded areas indicate LAT- 
detected (ELBS) quasars. 



source displays significant 7 -ray va riability during the first 3 
months of LAT observations (lAbdo et al. 2009b). In FigureH] 
we plot the distributions of our radio activity index Vq^ for 
three classes of sources, based on their LAT detection and 
LAT variability flag. The distributions indicate that radio- 
selected blazars were more likely to be detected in 7-rays by 
the LAT if they were in a particularly active radio state in May 
-December, 2008. 

The average radio activity index VJjg for the non-ELBS 
MOJAVE sources is -0.01 ±0.03. For ah ELBS-MOJAVE 
sources it is 0.23 ±0.07, and it rises to 0.31 ±0.09 if ELBS- 
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Fig. 4. — Distribution of the radio activity index V^aqj (see definition in 
the text) for the complete MOJAVE sample. The top histogram represents 
Fermi LAT non-detected (non-ELBS) sources within the sample. The bottom 
histogram represents the Fermi AGN 7-ray detections where black shading is 
used for objects with the 7-ray vaiiability flag from Abdo et al. ( 2009bj) . 
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MOJAVE sources with the LAT variability flag only are se- 
lected. The K-S test confirms, at a 99.0 per cent confi- 
dence, that the ELBS and non-ELBS MOJAVE sources have 
different distributions of VJfg. The confidence level rises to 
99.7 per cent if Vyg distribution for confirmed 7-ray variabile 
ELBS objects are compared to non-ELBS. Our conclusions 
stay the same even if a few outliers with very high radio ac- 
tivity index Vq^ > 1.0 are excluded from the analysis. In this 
case, the average Vq^ values decrease to 0.16 ±0.05 (ELBS- 
MOJAVE) and 0.24 ± 0.07 (ELBS-MOJAVE with the 7-ray 
variability flag). And we find the following confidence level 
of the K-S test for the VJg distribution: 97.2 per cent for ELBS 
vs. non-ELBS, 99.5 per cent for ELBS confirmed variable vs. 
non-ELBS. The fact that the average value of Vq^ for the entire 
MOJAVE sample (0.05 ± 0.03) is close to zero demonstrates 
a lack of overall bias. It is important to note that the same V 
activity analysis was performed for VLBA data prior to 2008 
and has shown no significant difference between ELBS and 
non-ELBS MOJAVE sources. 

Thus we observe a radio activity index significantly greater 
than zero for the ELBS sources only for the time interval of 
VLBA observations which roughly coincides with the ELBS 
detections. This temporal coincidence suggests that the 7- 
rays are produced in a region close to or coinc ident with the 
radio core where radio flares originate (e.g., iKovalev et al.l 
I2005h . However, more densely sampled 7-ray Fermi LAT — 
radio VLBI and/or single dish monitoring is necessary to es- 
tablish this robustly as the apparent time compression in rela- 
tivistic jets can be quite large. The frequency dependent opac- 
ity effect which delays radio flares should also be accounted 
for 

Our conclusion is supported by the fact that we find ex- 
tremely few blazars in their low radio activity state to be 7- 
ray bright. It is also interesting to note that the ELBS sources 
in Figure |4] with no variability flag (grey shading) span both 
moderate and high Vq^ states. This suggests that AGN in this 
category might consist of both true 7-ray non-variable objects 
as well as moderately flaring ones. 

3.3. Brightness temperature of jet cores 

Further insight into the 7-ray/radio jet connection can be 
obtained by examining the radio core brightness temperature 
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Fig. 5. — Distribution of median brightness temperature values, T\j of VLBI 
cores in tlie complete MOJAVE sample. Tlie sliaded areas represent tlie LAT- 
detected objects in the sample. 



7b. We have calculated median v alues of Tj, for the MOJA VE 
sample for the period 1999-2007 dMcCormick et al.ll2009l) . In 
less than 20 per cent of the cases only a lower limit could 
be found; this group of sources includes two ELBS objects. 
Gehan's generalized Wilcoxo n test fr om the ASURV sur- 
vival analysis package (Lavalley e t al.l fl 992) indicates, at a 
99.9 per cent confidence, that the median 7b values for LAT- 
detected sources are statistically higher than those for the rest 
of the sample (Figure|5]). The same result is found if the max- 
imum Ti, values are used in the analysis. Based on earlier 
brightness temperature findings of H oman et al. (2006), these 
results suggest higher Doppler fact ors for the LAT-detected 
sources, although their recent work dMcCormick et alj|2009l) 
indicates that variations in the intrinsic among jets may also 
play an important role in the brightest objects. Therefore, we 
expect that they might also have faster apparent jet speeds, 
which is indeed confirme d by results of the direct kinemat- 
ics analysis presented by iLister et al] (|2009b). Lister et al. 
have shown that LAT-detected quasars in the complete MO- 
JAVE sample have preferentially faster jet motions than the 
non-detected ones. 

4. SUMMARY 

On the basis of the joint analysis of the Fermi 7-ray LAT 
and radio observations of parsec-scale jets in blazars we con- 
clude the following. The 7-ray and parsec-scale radio emis- 
sions are strongly related in bright 7-ray objects detected by 
Fermi . At radio wavelengths, 7-ray bright sources are found 
to be preferentially brighter and more compact, which sug- 
gests that they might have higher Doppler factors than other 
blazars. The correlations found suggest that the prominent 
flares in both 7-ray and radio bands are produced in the cores 
of parsec-scale jets, typically within an apparent time separa- 
tion of up to a few months. These findings could be a conse- 
quence of relativistic beaming that boosts the jet emission in 
both bands. 

The first three months of Fermi observations represent a 
significant improvement over the earlier EGRET results, due 
to the dramatic increase in sensitivity and temporal coverage. 
The combination of these Fermi LAT measurements and the 
extensive VLBA monitoring at 15 GHz by the MOJAVE pro- 
gram has proved to be a powerful tool to study the nature of 
the emission processes in extragalactic jets. We found a clear 
connection between the beamed relativistic radio jets and the 
7-ray emission expected to originate in these regions, con- 
firming and enhancing earlier results obtained by comparing 
radio data with the EGRET catalog. 

We thank M. H. Cohen and A. R Lobanov for thorough 
reading of the manuscript and fruitful discussions. The au- 
thors wish to acknowledge the contributions of the rest of 
the MOJAVE team as well as students at Max-Planck-Institut 
fiir Radioastronomie and Purdue University. We also thank 
J. McEnery, D. Thompson and the Fermi LAT team for dis- 
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